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We demonstrate excitation of ferromagnetic resonance in CoFeB/MgO/CoFeB magnetic tunnel 
junctions (MTJs) by the combined action of voltage-controlled magnetic anisotropy (VCMA) and 
spin transfer torque (ST). Our measurements reveal that GHz-frequency VCMA torque and ST in 
low-resistance MTJs have similar magnitudes, and thus that both torques are equally important for 
understanding high-frequency voltage-driven magnetization dynamics in MTJs. As an example, we 
show that VCMA can increase the sensitivity of an MTJ-based microwave signal detector to the 
sensitivity level of semiconductor Schottky diodes. 
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Excitation of sub-nanosecond magnetic dynamics by 
electric field is a grand challenge in the field of spintron- 
ics. The ability to perform high-speed manipulation of 
magnetization by electric fields rather than by current- 
induced spin torques or magnetic fields would greatly re- 
duce ohmic losses and thereby improve the performance 
of spintronic devices such as non- volatile magnetic mem- 
ory. In this Letter we demonstrate excitation of GHz- 
range magnetization dynamics in nanoscale MTJs by the 
combined action of ST arising from spin-polarized cur- 
rent [1-9] and VCMA induced by electric field [10-20]. 
We show that ST and VCMA torques have similar mag- 
nitudes in MTJs with ultra-thin (< 1 nm) tunnel barriers 
for both DC and microwave-frequency voltages. Our re- 
sults demonstrate that a description of voltage-induced 
high-speed magnetic dynamics in MTJs should generally 
include not only ST but also VCMA terms. As an exam- 
ple, we show that VCMA can increase the sensitivity of a 
ST microwave signal detector, raising it to the sensitivity 
level of semiconductor Schottky diodes. 

We make measurements of voltage-driven magnetiza- 
tion dynamics in 150x70 nm 2 elliptical MTJ nanopil- 
lars shown schematically in Fig. 1. The nanopillars are 
patterned by ion milling from (bottom lead)/ Ta(5)/ 
PtMn(15)/ Co 70 Fe 30 (2.3)/ Ru(0.85)/ Co 4 oFe4oB 20 (2.4)/ 
MgO(0.83)/ Co 20 Fe 6 oB 2 o(1.58)/ Ta(5)/ (cap) multilayer 
(thicknesses in nm) with resistance-area product of 3.5 
f2 • /mi 2 deposited by magnetron sputtering in a Singu- 
lus TIMARIS system. Prior to patterning, the multilay- 
ers are annealed for 2 hours at 300 °C in a 1 Tesla in- 
plane magnetic field that sets the pinned layer exchange 
bias direction parallel to the long axis of the nanopillars 
(the ai-axis of our Cartesian coordinate system). The 
Co2oFegoB2o free layer exhibits perpendicular magnetic 
anisotropy [21], and the equilibrium direction of its mag- 
netization at zero field is normal to the sample plane 
(z-axis). Measurements of the nanopillar conductance G 
versus x-axis magnetic field H x (defined as the sum of 
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FIG. 1. (color online). Hysteresis loop of MTJ conductance 
G versus in-plane magnetic field H x measured under ± 80 
mV direct voltage bias Vd c - The left inset shows schematic of 
the MTJ nanopillar. The right inset displays variation of the 
effective perpendicular anisotropy field /\H p (Vd c ) with DC 
voltage extracted from the G(H X ) hysteresis loops via Eq.(2) 
(squares) and a linear fit to the data (line). 

the external .x-axis field and the stray field from the po- 
larizer) confirm the out-of-plane orientation of the free 
layer magnetization at H x = (see Fig. 1). The MTJ 
conductance exhibits cosine dependence on the angle <f> 
between the magnetic moments of the free and pinned 
layers [22]: 

G = G (1 + P 2 cos(^)), (1) 

where P is the current spin polarization. It follows that 
the shape of the G(H X ) curve is identical to that of the 
M X (H X ) = M s cos(0) hysteresis loop, where M s is satura- 
tion magnetization of the free layer and M x is projection 
of magnetization onto the x axis. The G(H X ) curves in 
Fig. 1 are hard-axis hysteresis loops, the slopes of which 
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FIG. 2. (color online). Rectified voltage V m ix versus fre- 
quency / of the applied microwave signal measured at sev- 
eral values of the in-plane magnetic field H x . The FMR re- 
sponse curves V m ix(f) at different fields are vertically offset by 
0.06 mV. At H x =0, Vmix(f) is a symmetric Lorentzian curve, 
while at \H X \ >0, V m ix{f) is a sum of symmetric and anti- 
symmetric Lorentzians. The antisymmetric part of V m i x (f) 
changes sign upon reversal of H x , while the symmetric part 
is negative for any H x . 



depend on the applied direct voltage bias Vd c by virtue 
of VCMA [17-20]. As discussed in the supplementary 
online material [23], ST from Vdc does not significantly 
affect the G(H X ) hysteresis loop. 

The effective perpendicular magnetic anisotropy en- 
ergy per unit volume of the free layer, E p , includes mag- 
netocrystalline, surface, and magnetostatic anisotropy 
contributions and can be calculated from the area un- 
der the hysteresis loop M X (H X ) in Fig. 1: 



1 



H x {M x )dM x = -M S H S 



(2) 



where H v is the effective perpendicular anisotropy field. 
The inset in Fig. 1 displays the variation of the perpen- 
dicular anisotropy field with DC voltage AH p (Vd c ) = 
H p (Vdc) — H p (0) . This dependence is well fit by a straight 
fine with the slope of 0.6 kOe/V, which corresponds to 
magnetic anisotropy energy per area per electric field of 
37 fJ/(Vm). This value is similar to the magnitude of 
VCMA reported by other groups [15, 17]. 

We study the effect of high-frequency voltage oscilla- 
tions on magnetization of the MTJ free layer using a 
ferromagnetic resonance (FMR) technique. In this tech- 
nique, magnetization oscillations are excited by a mi- 
crowave voltage applied to the MTJ, and the amplitude 
of magnetization precession is detected electrically via 
the rectified voltage V m i x generated by the MTJ [24, 25]. 
This technique is commonly called spin torque FMR 
(ST-FMR), since the magnetization precession of the 
free layer is driven by ST from the microwave current 



I ac cos(27r/i) flowing through the MTJ [26]. However, as 
we show below, magnetization oscillations in our MTJ 
samples are excited not only by ST from I ac but also 
by oscillating VCMA induced by the applied microwave 
voltage V ac - In this method, the rectified voltage V m i x = 
\l a cRac cos(V0 arises from mixing of the microwave cur- 
rent I ac cos(2tt ft) and tunneling magnetoresistance oscil- 
lations R ac cos(2tt ft + tp) resulting from the free layer's 
magnetization precession [25]. We sweep the microwave 
drive frequency / and measure ST-FMR response curve 
Vmi x (f). In general, peaks in the V m i x (f) response curve 
arise from resonant excitation of spin wave eigenmodes 
of the MTJ, and analysis of the line shape of these ST- 
FMR resonances gives the value of Gilbert damping and 
the magnitude of voltage-induced torques (including ST) 
acting on magnetization [26]. 

The measured dependence of V m i x (/) is shown in Fig. 2 
for several values of H x . The negative peaks in V m i x (f) 
arise from excitation of the quasi-uniform FMR mode of 
the free layer [25]. For all measurements, the applied 
microwave power is —36 dBm, which is small enough 
to excite the free layer precession in the linear regime. 
At H x = 0, the magnetization of the free layer is per- 
pendicular to the sample plane (0 = 0), and the ST- 
FMR line shape V m i x (f) is well described by a symmet- 
ric Lorentzian. For non-zero H x , the magnetization of 
the free layer is tilted away from the sample normal, and 
Vmix{f) develops asymmetry. In this regime, the ST- 
FMR line shape is well fit by a sum of symmetric and 
antisymmetric Lorentzians with identical resonance fre- 
quencies / : 
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where V s and V a are the symmetric and antisymmetric 
term amplitudes, and a is the half width at half maxi- 
mum. As is evident from Fig. 2 (compare ST-FMR curves 
at H x = ±120 Oe), the antisymmetric part of the ST- 
FMR line shape, V a , changes sign upon reversal of H x . 

The symmetric component of the ST-FMR curve arises 
from in-plane ST m x (e p x rh) that lies in the plane 
defined by the magnetic moments of the free and the 
pinned layers, where m is the unit vector in the direc- 
tion of the free layer's magnetization and e p is the unit 
vector in the direction of the polarizer's magnetization. 
An antisymmetric component of the ST-FMR curve has 
been previously observed in nanoscale MTJs at non-zero 
DC bias voltages, and was shown to arise from a field-like 
ST, Ty~ rh x e p [27, 28]. The antisymmetric component 
was observed to vanish at zero bias voltage, which was 
explained by quadratic dependence of Tf on the voltage 
bias [27]. Hence, our observation of a non-zero antisym- 
metric component of the ST-FMR response curve at zero 
DC voltage is surprising. 

In order to understand the origin of the ST-FMR line 
shape asymmetry in our MTJ system, we numerically 
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FIG. 3. (color online). Macrospin simulations of ST-FMR 
line shape assuming non-zero in-plane ST and (a) non-zero 
VCMA (A22± = 4 Oe) and zero field-like ST (b f = 0), (b) 
zero VCMA (AH x =0 Oe) and non-zero field-like ST (&/= 
0.5). The curves for different values of H x are vertically offset 
by 0.06 mV. 



simulate [23] the resonance curve by solving the Landau- 
Lifshitz-Gilbert equation of motion of the free layer in 
the macrospin approximation [29]: 

dm „ „ dm. Jh 

— = - 7 mxff eff + amx— - 7 ^— 



1 + P 2 cos(</>) 



[m x (m x e p ) + bjrh x e p ]. (4) 



In Eq.(4), H e ff is the effective field including contribu- 
tions from the external, magnetostatic, and perpendic- 
ular magnetic anisotropy fields, 7 is the gyromagnetic 
ratio, a is the Gilbert damping parameter, Jh/2e is the 
microwave spin current density, d is the thickness of the 
free layer, and bf is the ratio of the magnitudes of the 
in-plane ST and field-like ST. 

For the approximately quadratic dependence of 77 
on the current bias observed in previous experiments 
[27, 28, 30], bf is proportional to J and thus 77 can be 
neglected at zero current bias. However, we include 77 in 
the simulations since a linear dependence of 77 on current 
(corresponding to bf = const) is, in principle, possible for 
MTJs with different free and pinned ferromagnetic lay- 
ers [31]. Such a linear dependence of 77 on current would 
give rise to an antisymmetric component of the ST-FMR 
response curve at zero bias. 

In our simulations, both the charge current density 
J = Jo sin(27r/i) and the effective magnetic field f/ e ff 
are generally time-dependent. The time-dependence of 
the effective magnetic field arises from the VCMA, giv- 
ing the time-dependent perpendicular anisotropy field 
H ±z = (H ±0 + AH ± sin(27r/i)) cos(0), where H ±0 is the 



anisotropy field at zero bias and AH± is the VCMA 
field. The simulated ST-FMR response curves shown 
in Fig. 3(a) are calculated assuming non-zero in-plane 
ST and non-zero VCMA (AH ± = 4 Oe), but zero field- 
like ST (bf = 0). In these simulations, we use H± Q = 
11,362 Oe, a = 0.030, and P = 0.54 as determined by 
fitting the resonance frequency, the line width and the 
amplitude of the ST-FMR resonance curve measured at 
H x = 0. Allowing for uncertainty in the sample vol- 
ume, we note that P = 0.54 agrees well with the value 
P = 0.52 calculated from the MTJ conductance using 
Eq.(l). In the simulations, we use demagnetization fac- 
tors N x = 0.014, N y = 0.040, N z = 0.946 obtained from 
the assumption that the free layer is an elliptic cylin- 
der [32] . The saturation magnetization M s is taken to be 
950 emu/cm 3 as measured by vibrating sample magne- 
tometry, and AH± = 4 Oe is calculated from the slope 
of AH p (Vdc) line in the inset of Fig. 1. Fig. 3(a) shows 
that the simulated ST-FMR line shapes reproduce the 
symmetry of the measured curves: (i) the line shape is 
symmetric at H x = and (ii) the line shape is asym- 
metric with the antisymmetric part changing sign upon 
reversal of H x for \H X \ > 0. In contrast, simulations with 
a non-zero filed-like torque 77 do not reproduce the sym- 
metry of the measured line shapes. Fig. 3(b) shows the 
simulation results for non-zero in-plane ST and non-zero 
field-like ST {b f = 0.5) but zero VCMA (AH ± = 0). The 
simulated ST-FMR line shapes are asymmetric for all ap- 
plied field values, and the sign of the antisymmetric part 
does not depend on the applied field direction. Therefore, 
our simulations demonstrate that the observed ST-FMR 
line shape asymmetry arises from VCMA, not field-like 
ST. 

The antisymmetric part of the ST-FMR line shape cal- 
culated with the VCMA value derived from DC measure- 
ments (AHj_ = 4 Oe) has similar magnitude to that of 
the measured ST-FMR response curve. This implies that 
the magnitude of VCMA is comparable for microwave 
and DC voltages. Since the measured magnitudes of the 
symmetric and antisymmetric parts of the ST-FMR line 
shapes are similar for 0^0, we conclude that ST and 
VCMA torque generally have similar magnitudes in our 
MTJs. 

Although the ST-FMR line shape is well reproduced 
by the simulations with the VCMA drive term, the de- 
pendence of the simulated resonance frequency /n on H x 
does not qualitatively agree with the measurements. The 
macrospin simulations predict decreasing /o with increas- 
ing 1 22x1 as shown in Fig. 4 (a,c), while the measured 
fo{\Hx\) shows the opposite trend (Fig. 2). We find that 
this discrepancy arises from the simplified form of the 
perpendicular anisotropy energy used in the simulations. 
The thickness of the free layer in our samples is intention- 
ally tuned to a value at which the uniaxial perpendicu- 
lar anisotropy and the easy-plane shape anisotropy have 
comparable magnitudes but opposite signs. Therefore, 
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FIG. 4. (color online). Macrospin simulations of ST-FMR 
response curves with (a) zero (-ffxi = 0) and (b) positive 
(H±\ = 280 Oe) second-order perpendicular anisotropy. The 
curves are vertically offset by 0.06 mV. (c) Macrospin and mi- 
cromagnetic simulations of ST-FMR resonance frequency ver- 
sus H x for zero and positive Hj_i (lines) shown together with 
the experimental data (solid squares), (d) Macrospin and mi- 
cromagnetic simulations (lines) of the G(H X ) hysteresis loop 
using the magnetic anisotropy parameters given by the best 
fits in (c). Squares are experimental data measured at zero 
bias voltage, (e) Sensitivity of the MTJ to a microwave signal 
is enhanced by VCMA when the free layer magnetization is 
tilted from the sample normal by magnetic field H x = +40 
Oe. 

the leading terms proportional to sin 2 (f9) in the angular 
dependence of both anisotropy energy expansions nearly 
cancel each other and the higher order term proportional 
to sin 4 (f?) becomes important [33]. 

In order to study the impact of the second-order 
perpendicular anisotropy term on fo(H x ), we repeat 
our macrospin simulations with non-zero first- and 



second-order perpendicular anisotropy terms: E± = 
K\ sin 2 (0) + I\ 2 sin 4 (6>), which gives H± z = [Hxo + 
H±i(l — cos 2 (f9))] cos(6>) at zero bias, where H±q = 
2K\jM s and H±i = AK 2 /M S . Our simulations show 
that adding a small positive second-order anisotropy field 
H±i changes the shift of the resonance frequency with 
\H X \ from negative to positive, and the results from sim- 
ulations at H±i = 280 Oe are in good quantitative agree- 
ment with the measured fo(H x ) (Fig. 4 (b,c)). 

We also make micromagnetic simulations of fo(H x ) us- 
ing MicroMagus simulator and find the micromagnetic 
fo(H x ) to be in qualitative agreement with the macrospin 
simulation results. Fig. 4 (c) shows that micromagnetic 
simulations performed with H±i = give fo(H x ) in- 
consistent with the experimental data, while those per- 
formed with Hxo = 11,515 Oe and H ±1 = 240 Oe re- 
produce the measured fo(H x ) well. Micromagnetic and 
macrospin simulations of the G(H X ) hysteresis loop made 
with the values of H± and H±i obtained from best fits of 
fo(H x ) are also in good agreement with the experimental 
data as shown in Fig. 4(d). 

The measured amplitudes of the symmetric and anti- 
symmetric parts of the ST-FMR line shape are similar 
to each other for \H X \ > 0. This implies that the mag- 
nitude of VCMA torque is similar to that of ST, and 
thus both torques are important for the description of 
voltage-driven magnetization dynamics in MTJs, such 
as ST- induced magnetization reversal [34, 35]. As an 
example, Fig. 4 (e) shows that VCMA can significantly 
increase the sensitivity of an MTJ-based microwave sig- 
nal detector [24]. Indeed, the maximum rectified voltage 
generated by the MTJ in response to applied microwave 
signal increases by 39% when H x increases from Oe 
to +40 Oe [23]. This increase in the detector sensitiv- 
ity is mainly due to the antisymmetric part of the ST- 
FMR line shape, which is induced by VCMA. With the 
VCMA contribution, the MTJ microwave detector sen- 
sitivity reaches 440 V/W — a value similar to that of 
semiconductor Schottky diodes (~ 5 x 10 2 V/W [36]). 

In conclusion, we demonstrate excitation of GHz-range 
magnetization dynamics in nanoscale MTJs by the com- 
bined action of spin torque and voltage-controlled mag- 
netic anisotropy. Our work shows that the magnitudes of 
high-frequency spin torque and voltage-controlled mag- 
netic anisotropy torque in MgO-based MTJs can be sim- 
ilar to each other, and thus that quantitative descrip- 
tions of voltage-driven magnetization dynamics in MTJs 
should generally include both torque terms. We show 
that voltage-controlled magnetic anisotropy can increase 
the sensitivity of spin torque MTJ microwave signal de- 
tectors to the sensitivity level of Schottky diode detec- 
tors. 
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Supplementary Material 

1. Macrospin Simulations 

We simulate dynamics of the MTJ free layer by numerically solving Landau-Lifshitz-Gilbert (LLG) 
equation with the spin torque term in the macrospin approximation: 

dm „ - „ dm Jh P ,„„.,„„ n .... 

= -ym x H /lff +amx y [mx(mxe„) + b f mxe„], (S1) 

dt 11 dt / 2eM s d l + P 2 cos(0) p f p 

where m{t)\s the unit vector in the direction of the free layer magnetic moment. In the LLG equation, 

j(t) = J cos(2^?), where /is the frequency of the applied microwave voltage, H eff = H exl +H Jemag +H p \s 

the total effective field including the external H ext =H x x , demagnetizing 

H demag =-N x AnM s m x x-N y AnM s m y y-N 2 AnM s m 2 z and perpendicular magnetic anisotropy 

H = (H ± + AH L sm(2nft))m z z fields. Here Hjs the static perpendicular anisotropy field at zero voltage 

bias and AH ± sin(2^) is the time-dependent perpendicular anisotropy field arising from the applied ac 
voltage via VCMA. The value of AH ± = 4 Oe is derived from the change of slope of the hard axis 
hysteresis loop with applied dc voltage as shown in Fig. 1 . The angular dependence of the perpendicular 



anisotropy field H ± m z =[H 10 +H n (l-cos 2 6>)]cos& is derived from the angular dependence of the 

perpendicular magnetic anisotropy energy density^ , in which first- and second-order terms are included: 

sin 4 9 , where H L0 = 2K l I M s and H L1 = 4K 2 1 M s ■ The magnitude of saturation 

magnetization of the free layer M s = 950 emu/cm 3 is directly measured by vibrating sample magnetometry 
of unpatterned multilayer annealed under the same conditions as the multilayer used for the nanopillar 
fabrication. We employ the demagnetization factors N x = 0.014, N y = 0.040, N z = 0.946 appropriate for the 
150 x 70 x 1 .58 nm 3 elliptical cylinder, which are the dimensions of the free layer. The value of the Gilbert 
damping constant a = 0.030 is determined by fitting the line width of the measured ferromagnetic 
resonance curve while the degree of spin polarization of the current P = 0.54 is determined from fitting the 
amplitude of the ferromagnetic resonance curve. The magnitude of the first-order perpendicular 
anisotropy field H ±0 = 1 1 ,362 Oe is given by the best fit of the resonance frequency at zero field (H x = 0). 

The magnitude of the second-order perpendicular anisotropy field H u = 280 Oe is given by the best fit of 
the dependence of the resonance frequency on applied magnetic field frfHJ. 

We numerically solve the LLG equation for a period of time sufficiently long for the free layer 
magnetic moment m{t) to reach the regime of steady-state oscillations. In this steady state regime, we 

use m(t) to calculate the time-dependent MTJ conductance G = G (l + P 2 cos(^))= G (l + P 2 m x (t)), and tne 

time-dependent voltage generated by the MTJ due to the conductance oscillations^ (A = ^9 cos ( 2 #) . 

* W G {l + P 2 m x (t)} 

The mixing voltage generated by the MTJ, V mix , at the drive frequency /is calculated as time average of 
V g (t) over one period of the magnetization precession. Performing the calculation of V mix as a function of 
the drive frequency f, we obtain the simulated ferromagnetic resonance curve V mix (f) in the macrospin 
approximation. 

2. Effect of direct voltage on the hysteresis loop of the free layer 

Direct voltage applied to the MTJ generally induces in-plane spin torque, field-like torque and 
VCMA torque, each of which can potentially alter the hysteresis loop of magnetization versus magnetic 
field M/HJ. However, each type of torque is expected to produce a distinctly different effect on the M X (H X ) 



hysteresis loop and thus measurements of M X (HJ as a function of voltage can give us information on the 
type of voltage-induced torques present in the system. 

Effective magnetic field arising from the in-plane spin torque is perpendicular to the plane defined 
by the magnetic moments of the free layer and the polarizer. At zero applied voltage, magnetic moment of 
the polarizer is parallel to the x-axis and the free layer is in the x-z plane, where z is the sample normal. 
Therefore, the effective field form the in-plane spin torque is in the >> direction, which gives rise to rotation 
of the free layer magnetization in the direction perpendicular to the x-z plane. As a result, the angle 
between the magnetic moments of the free layer and the polarizer <j> (and thus M x = M s cos(<j>)) remains 
unchanged to first order in the free layer rotation angle. For this reason, we do not expect the in-plane 
spin torque to significantly modify the M X (H X ) hysteresis loop. 

The effective field arising from the field-like torque is collinear with the polarizer direction (along 
the x-axis), which is the same as the direction of the applied magnetic field H x . As a result, the field-like 
torque is expected to shift the hysteresis loop along the magnetic field axis but it cannot change the hard- 
axis hysteresis loop slope. Since we do not see a significant voltage-induced shift of the measured 
hysteresis loop along H x , we conclude that field-like torque is relatively small for the range of direct 
voltages employed in the experiment. 

The effective field arising from VCMA is parallel to the z-axis and it causes the magnetization of 
the free layer to rotate in the x-z plane thereby changing the angle <j> between the free layer and the 
polarizer (and thus M x = M s cos(</>)). When we sweep the external magnetic field, the polar angle of the free 
layer will change and hence change the effective field. Therefore, the shift of the hysteresis loop is not a 
constant as opposite to the field-like torque. In fact, it will change the slope of the hysteresis loop. That 
explains that the change of slopes of the hysteresis loop is due to the VCMA. 

3. Field-dependence of the MTJ microwave detector sensitivity 

The sensitivity of a MTJ-based microwave signal detector is defined asK A . I P tn , where V dc is 
direct voltage generated by the MTJ and P in is the applied microwave power. In our experiment, the 
rectified voltage V mjx = (y[2 1 n)V dc is measured with a lock-in amplifier. Fig.SI shows the measured 



dependence of the MTJ sensitivity on in-plane magnetic field H x . The maximum sensitivity is observed at 
H x > due to a non-zero VCMA torque acting on the free layer magnetization tilted away from the sample 
normal by H x . 




FIG.S1 . Sensitivity of the MTJ as a function of in-plane magnetic field. 



